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The A cold dark matter (ACDM) model is currently known as the simplest cosmology model that 
best describes observations with minimal number of parameters. Here we introduce a cosmology 
model that is preferred over the conventional ACDM one by constructing dark energy as the sum 
of the cosmological constant A and the additional fluid that is designed to have an extremely short 
transient spike in energy density during the radiation-matter equality era and the early scaling 
behavior with radiation and matter densities. The density parameter of the additional fluid is defined 
as a Gaussian function plus a constant in logarithmic scale-factor space. Searching for the best-£t 
cosmological parameters in the presence of such a dark energy spike gives a far smaller chi-square 
value by about hve times the number of additional parameters introduced and narrower constraints 
on matter density and Hubble constant compared with the best-fit ACDM model. The signihcant 
improvement in reducing chi-square mainly comes from the better fitting of Planck temperature 
power spectrum around the third {£ « 800) and sixth {£ « 1800) acoustic peaks. The likelihood ratio 
test and the Akaike information criterion suggest that the model of dark energy spike is strongly 
favored by the current cosmological observations over the conventional ACDM model. However, 
based on the Bayesian information criterion which penalizes models with more parameters, the 
strong evidence supporting the presence of dark energy spike disappears. Our result emphasizes that 
the alternative cosmological parameter estimation with even better fitting of the same observational 
data is allowed in the Einstein’s gravity. 

PACS numbers: 98.80.-k, 95.36.+x 


I. INTRODUCTION 

It is one of the primary aims in cosmology to find the 
simplest model that best describes the astronomical ob¬ 
servations. Until now, the best concordance model with 
minimal number of parameters is currently known as the 
A cold dark matter (ACDM) model with the cosmolog¬ 
ical constant A as dark energy and CDM as dominant 
dark matter [in. 

Many efforts have been made to develop a cosmology 
model that is better than the conventional ACDM model. 
Determining which model is preferred over the other is 
a problem of model selection. The standard approach 
used in the model selection is that one constrains the 
new model with data using the likelihood method and 
checks whether or not it is supported over the previous 
best model based on the statistical criteria. In usual 
cases, the new candidate cosmology model has more free 
parameters than ACDM model while it often gives a bet¬ 
ter fitting of observational data. However, simply adding 
more parameters and getting smaller chi-square (or larger 
likelihood) does not make the relevant model rank as a 
better model. In order for a new cosmology model to be 
ranked as a model better than ACDM model, it should 
pass through at least one of the model selection criteria 
such as the likelihood ratio test Q, Akaike information 
criterion Q , Bayesian information criterion Q , Bayesian 
evidence, and so on (see also [Ml for applications in 
cosmology with brief reviews). 

Although many dark energy models have been pro¬ 
posed, most of them give only a small improvement in 
fitting the observational data compared with the ACDM 


model and do not pass through the model selection cri¬ 
teria with high significance, implying that the ACDM 
model is the final winner in the competition of model 
selection (Tol - [^ . 

Recently, Park et al. [11 investigated the observational 
effect of the early episodically dominating dark energy 
based on the minimally coupled scalar field with the 
Albrecht-Skordis potential, where the dark energy den¬ 
sity transiently becomes strong during a short period of 
time. They show that the presence of the early epdisodic 
dark energy can affect the cosmological parameter esti¬ 
mation significantly and conclude that the recent Planck 
data strongly favor the ACDM model because only a 
limited amount of dark energy with episodic nature is 
allowed. In this paper, we introduce a fluid version of 
early transiently dominating dark energy model with the 
similar episodic nature. Our dark energy model is de¬ 
signed to have a transient spike in energy density during 
an extremely short period and the early scaling behavior 
with radiation and matter density. We show that our 
dark energy model gives a significant improved fit to the 
recent observational data with different parameter con¬ 
straints and thus is preferred over the best-ht ACDM 
model based on some model selection criteria. Through 
the example of dark energy spike model, we show that the 
alternative parameter estimation with even better fitting 
of the same observational data is allowed in the Einstein’s 
gravity. 

This paper is organized as follows. Sec. ITT] describes 
the fluid-based dark energy spike model with a transient 
variation in dark energy density and presents numerical 
calculations of background evolution of this model. In 
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Sec. mi observational effects of the dark energy spike 
are investigated using the recent observational data such 
as the cosmic microwave background radiation (CMB) 
data from the Planck satellite and the baryonic acous¬ 
tic oscillation (BAO) data from the large-scale structure 
surveys. The cosmological parameters constrained with 
observations are compared in the presence of or with¬ 
out the dark energy spike. In Sec. lIVi we compare our 
dark energy model with the conventional ACDM model 
based on some statistical criteria used in model selec¬ 
tion. The discussion and conclusion are presented in Sec. 
rvl Throughout this paper, we set c = 1 and SttG = 1. 


II. A FLUID-BASED DARK ENERGY SPIKE 
MODEL 


The quintessence-based early episodically dominating 
dark energy model proposed in [22j| is on a solid theoreti¬ 
cal footing, but has its limitations. First, it is not easy to 
control the onset, strength, and duration of the transient 
dark energy because the behavior of dark energy strongly 
depends on potential parameters and initial conditions. 
Second, the scalar-field based dark energy model theo¬ 
retically does not accommodate the crossing of phantom 
divide {w = —1) in the dark energy equation of state, 
and thus it is not allowed to consider a spike, a transient 
and abrupt variation of dark energy density, which in¬ 
evitably induces re = — 1 crossing. Here we introduce a 
fluid model of early dark energy that allows a dark energy 
spike and is easy to handle numerically. 

In the conventional ACDM model, we add a new fluid 
(denoted as x) with a transient spike in energy density. 
We assume that the dark energy density parameter (VLx) 
is represented as the sum of a Gaussian function and a 
constant in logarithmic scale-factor space. 


^x<,a) = = Aexp 


(In a — In 
2^2 


+ B, 


( 1 ) 


the ACDM model (A = 0, i? = 0) and the scaling-DE 
model (A = 0 and B > 0) are nested within the spike-DE 
model. 

In the presence of the x-fluid, the squared Hubble pa¬ 
rameter normalized with the present value is given by 

f B \ _ DrOO ^ -|- ^mOO. ^ -|- Dao + ^KQO. ^ 

\Jh) ~ 1-D,(a) ’ ^ ^ 


where the subindices r, m, K represent the radia¬ 
tion, matter, and spatial curvature, respectively, and 
the subindex zero indicates the present value. We as¬ 
sume that the x-fluid satisfies the continuity equation, 
fix = —3H{p,x +Px)- The pressure of the x-fluid is given 

by 


Px = 3H^nx -1 - 
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Throughout this paper, we consider the spatially flat uni¬ 
verse (Da-o = 0). The equation of state of the effective 
dark energy fluid (denoted with a subindex d) becomes 


w = 


P£ 


2HVLx 
+ ^a) 


^'x 

3(172, -|- Da) 


( 6 ) 


where fix is the energy density of the additional fluid, a{t) 
is the cosmic scale factor normalized to unity at present, 
H = d/a is the Hubble parameter at epoch a with a dot 
as a time-derivative, and A, Oc, and a are interpreted as 
amplitude (strength), temporal position, and duration of 
the dark energy spike, respectively. The constant term B 
denotes the level of early dark energy that exists from the 
beginning of the universe and is responsible for the scal¬ 
ing evolution where the dark energy density follows that 
of the dominant fluid. In this paper, we incorporate the 
x-fluid and the cosmological constant A into the effective 
fluid of dark energy (DE). The behavior of dark energy 
in our model shows early scaling evolution with radiation 
and matter densities with a sudden dark energy spike at 
a particular epoch and the late-time acceleration phase 
due to the cosmological constant. From here on, we call 
the dark energy model with both scaling and transient 
behaviors as spike-DE model, and the model only with 
the scaling behavior as scaling-DE model. Note that both 


where fid = fJ-x + A, Pd = Px - A, and Da = 

The dark energy density parameter (D^j = Dj, -I- Da) has 
three asymptotic values: D^ ~ A -|- H at the onset of the 
dark energy spike (a = Oc), fid — B before and after the 
onset (|a —Ocl ^0) and before the late-time acceleration 
(a oq), and D^ ~ H-I-Da during the acceleration phase 
(a w ao > Oc). 

Eigure[T]shows the evolution of density parameters, en¬ 
ergy densities, dark energy equation of state, and Hubble 
parameter in the spike-DE model where a strong dark 
energy spike occurs at a = 10“^’^ (InUc = —8.059), to¬ 
gether with those in scaling-DE and ACDM models. As 
expected, the dark energy densities of spike-DE (green) 
and scaling-DE (grey) models show scaling behaviors fol¬ 
lowing radiation and matter sequentially. The model 
parameters have been adopted as the best-fit ones ob¬ 
tained with the recent observational data (see Sec. IHII 
for details). In the spike-DE model, the dark energy 
equation of state experiences a change of about three 
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FIG. 1: Evolution of density parameters (fli; i = r, m, d), energy densities (/ri), dark energy equation-of-state {w), and 
normalized Hubble parameter divided by the ACDM prediction [{H/Ho)/{H/Ho)acdm] in the best-fit spike-DE model with 
logj^Q ac = —3.5, A = 0.28, B = 1.1 x 10”"^, a = 1.5 x 10~^. In the top-panels, the behaviors of radiation (r), matter (m), and 
dark energy (d) components are shown as red, yellow, and green curves, respectively. The energy density due to the cosmological 
constant in the best-fit ACDM model is shown as a blue curve (top-right panel). In all panels, grey curves represent the results 
of scaling-DE model with B = 5.2 x 10“®. The precise values of model parameters used in the numerical calculation are 
presented in Table |T] 


orders of magnitude with the crossing of phantom di¬ 
vide twice during the occurrence of dark energy spike. 
Considering 95.4% (2cr) conhdence limits of the Gaus¬ 
sian shape of the spike, it lasts for about 600 years 
(3.153 X 10-4 <a< 3.172 x IQ-^). 


III. OBSERVATIONAL CONSTRAINTS ON THE 
DARK ENERGY SPIKE MODEL 

We probe the observational signatures of our spike-DE 
model by considering both the scalar- and tensor-type 
perturbations in a system of multiple components for ra¬ 
diation, matter, and effective dark energy fluid (the x- 
fluid plus the cosmological constant). For this aim, we 
have modified the publicly available CAMB/CosmoMC 
package (version of Dec. 13 2013) [l^, to include the 
evolution of background and perturbation of the effec¬ 


tive dark energy fluid, and explored the allowed ranges 
of the conventional cosmological parameters in the pres¬ 
ence of a dark energy spike using the Planck CMB data 
together with other external data sets. For the evolu¬ 
tion of perturbed density and velocity of the x-fluid, the 
parametrized post-Friedmann prescription for the dark 
energy perturbations is used to allow the multiple cross¬ 
ing of phantom divide {w = —1) in the time-dependent 
dark energy equation of state . Following the Planck 
team’s analysis, we assume the current CMB tempera¬ 
ture as To = 2.7255 K and the effective number of neu¬ 
trinos as = 3.046 with a single massive eigenstate of 
mass nil, = 0.06 eV 0. 

We use the CMB data obtained with Planck 37| , which 
is a combination of the CMB temperature anisotropy an¬ 
gular power spectrum up to small angular scales (£ = 
2500) and the Wilkinson Microwave Anisotropy Probe 9- 
year polarization data [2^. We used four Planck CMB 
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TABLE I: Best-fit cosmological parameters of the spatially 
flat ACDM, scaling-DE, and spike-DE models. 


Parameter 

ACDM 

Scaling-DE 

Spike-DE 

A 

0 

0 

0.28360 

B 

0 

5.1736 X 10“® 

1.0662 X 10“'‘ 

logic 



-3.5 

a 



1.4604 X 10“® 


2.21632 

2.22226 

2.19523 

Q.ch? 

0.11827 

0.11778 

0.11777 

h 

0.67929 

0.68130 

0.68158 

T 

0.09623 

0.08908 

0.08873 

Us 

0.96550 

0.96487 

0.97147 

r 

0.00048 

0.00046 

0.01640 

ln[10i°A,] 

3.09985 

3.08397 

3.08691 

to (Gyr) 

13.7992 

13.7930 

13.7983 

Aioo 

178.3636 

138.0731 

140.2106 

Aria 

62.92783 

50.31658 

61.88821 


118.6188 

115.4876 

126.9649 

aCIB 

^143 

6.620212 

3.852115 

5.640900 

aCIB 

^217 

25.52911 

26.94230 

23.30611 

AtSZ 

^143 

3.724382 

8.408760 

2.995350 

„PS 

^143x217 

0.9075909 

0.8956619 

0.9206412 

«CIB 

'143x217 

0.2190109 

0.3866272 

9.2171337 X 10”^ 

,yCIB 

0.5448702 

0.5265283 

0.5609424 

ClOO 

1.000590 

1.000599 

1.000575 

C217 

0.9963431 

0.9962796 

0.9968735 

^tSZ-CIB 

0.5315524 

5.734012 X lO"'^ 

0.2737207 

^kSZ 

0.1122116 

0.2684244 

0.1527581 

/3i 

0.5376251 

0.5772729 

0.2189442 


likelihood data sets (version of 2013), Lowlike for low £ 
temperature and polarization likelihood covering £ = 2- 
32, Commander for low £ temperature-only likelihood cov¬ 
ering £ = 2-49, CamSpec for high £ temperature-only like¬ 
lihood with £ = 50-2500 [l^, and Leasing for lensing 
effect [ 2 ^. As the external data derived from the large- 
scale structure observations, we also have used the BAO 
data points measured by Sloan Digital Sky Survey Data 
Release 7 (DR7) [ 2 ^, Baryon Oscillation Spectroscopic 
Survey Data Release 9 (DR9) 6dF Galaxy Survey 
[lOl, and Wiggle Z surveys [32 |. 

With CMB and BAO data, we have constrained the 
parameter space of the spatially flat ACDM, scaling-DE, 
and spike-DE models that are favored by the observa¬ 
tions. We limit our investigation by considering a spike- 
DE model with a dark energy spike occurring near the 
radiation-matter equality era (logj^g Oc = —3.5); see Fig. 
[TJ The reason for choosing such an epoch is that the 
transient domination of dark energy near the radiation- 


TABLE II: Mean and standard deviation (68.3% confidence 
limit) of the conventional cosmological parameters estimated 
from the marginalized one-dimensional likelihood distribution 
for best-fit ACDM, scaling-DE, spike-DE models constrained 
with the Planck CMB and BAO data sets. For the tensor-to- 
scalar ratio r and the level of early dark energy B, the upper 
limits are presented. 


Parameter 

ACDM 

Scaling-DE 

Spike-DE 

100^^60 

4.784 ± 0.076 

4.870 ±0.102 

4.727 ± 0.026 

Dco 

0.2547 ±0.0081 

0.2564 ± 0.0084 

0.2518 ±0.0031 

h 

0.6807 ± 0.0069 

0.6761 ± 0.0077 

0.6836 ± 0.0037 

T 

0.090 ± 0.012 

0.093 ±0.013 

0.086 ±0.011 

Us 

0.9644 ± 0.0056 

0.9675 ± 0.0061 

0.9635 ± 0.0046 

r 

< 0.054 

< 0.061 

< 0.048 

ln[10i“A.] 

3.084 ± 0.022 

3.088 ± 0.024 

3.081 ±0.021 

to (Gyr) 

13.795 ± 0.036 

13.858 ± 0.060 

13.781 ± 0.033 

B 


< 0.0022 



matter equality most affects the evolution of density 
perturbations in the scalar-field-based early dark energy 
model, inducing a highly oscillatory feature in the angu¬ 
lar power spectrum of temperature fluctuations at high 
multipoles (see Fig. 3 of Ref. [1^). The free conventional 
cosmological parameters are h, t, n^, r, 

and ln[10^°As], where Dbo (Idco) is the baryon (CDM) 
density parameter at the current epoch, h is the normal¬ 
ized Hubble constant with Hq = lOOh km s“^ Mpc“^, t 
is the reionization optical depth, Us is the spectral index 
of the primordial scalar-type perturbation, r is the ra¬ 
tio of tensor- to scalar-type perturbations, and As is the 
amplitude of the primordial curvature perturbations with 
As = Pfi{k)/(27:^) at the pivot scale fcg = 0.05 Mpc“^. 
The running of spectral index is not considered. There 
are also several foreground and calibration parameters 
(see Table [T] and 0 for detailed descriptions of the pa¬ 
rameters) . 

The free parameters of spike-DE model are A, B, and 
tr with Oc hxed [see Eq. ©]■ With conventional and dark 
energy model parameters all freely varying, the Markov 
chain Monte Carlo (MCMC) chains are not easily con¬ 
verged due to multiple local maxima in the multidimen¬ 
sional likelihood distribution. In this work, instead of ob¬ 
taining the full converged MCMC chains for all the free 
parameters, we search for the best-ht location in the like¬ 
lihood distribution by manually running the CosmoMC 
with an option action=2 starting at the local maxima 
found from the trial MCMC chains obtained with all pa¬ 
rameters varying. The results are summarized in Ta¬ 
ble [H which lists the parameters of ACDM, scaling-DE, 
and spike-DE models that best describe the observational 
data, together with the cosmic age (£o) and the param¬ 
eters related with foregrounds and instrumental calibra¬ 
tions. 

To see how the conventional cosmological parameters 
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FIG. 2: Top: Two-dimensional likelihood contours of the conventional cosmological parameters favored by the Planck CMB and 
BAO data sets for spike-DE (red), scaling-DE (grey), and ACDM (black curves) models. For the spike-DE model, dark energy 
parameters have been fixed with the values given in Table [H The thick and thin solid curves indicate the 68.3% and 95.4% 
confidence limits, respectively. Middle and bottom: Marginalized one-dimensional likelihood distributions for each cosmological 
parameter, with arbitrary normalizations. 


are affected by the presence of the dark energy spike, 
we apply the MCMC method to randomly explore the 
parameter space that is favored by observations. For 
the spike-DE model, we have fixed dark energy param¬ 
eters A, B, and a with the best-fit values given in Ta¬ 
ble ID For the scaling-DE model, however, the parameter 
B, the initial level of early dark energy, has been freely 
varied. Figure [5] shows two-dimensional likelihood con¬ 
tours and marginalized one-dimensional likelihood distri¬ 
butions of conventional cosmological parameters favored 
by the Planck CMB and BAO data sets for the spike- 
DE, scaling-DE, and ACDM models, estimated from the 
converged MCMC chains. Note that here we present the 
likelihood distributions of Q,bo and Dco instead of flboh'^ 
and Table llll summarizes the mean and 68.3% 

confidence limit of cosmological parameters estimated 
from the one-dimensional likelihood distributions. For 
the tensor-to-scalar ratio r and the level of early dark en¬ 
ergy B, the upper limits (68.3%) are given. Interestingly, 
compared with ACDM model, the spike-DE model gives 
narrower parameter constraints on baryon, CDM den¬ 
sity parameters and Hubble constant with the standard 
deviations smaller by a factor of 2.9, 2.6, 1.9, respec¬ 


tively, and best-fit values slightly deviating from those 
of ACDM model. Since the likelihoods of the spike-DE 
model sufficiently overlap with the ACDM ones, the es¬ 
timated parameters of both models are still consistent 
with each other. 

For the scaling-DE model, the parameter constraints 
are consistent with those of ACDM model, except for 
slightly larger values of baryon density and cosmic age. In 
this model we have set B as a free parameter to constrain 
the level of early dark energy density (He = B). The al¬ 
lowed range for the early dark energy is He < 0.0045 
(95.4% confidence limit), which is narrower than the 
Planck constraint on the fluid-based early dark energy 
density parameter of Doran & Robbers Q (De < 0.009; 
i)- Recently, a substantial improvement on the con¬ 
straint De < 0.0036 (at 95% confidence level) has been 
obtained by the Planck 2015 data analysis [^ . 

IV. MODEL COMPARISON 

In this section, we compare the best-fit ACDM, scaling- 
DE, spike-DE models to see which model is preferred 
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TABLE III: Chi-square (x^) values of best-fit ACDM, scaling- 
DE, and spike-DE models, together with differences of chi- 
square (Ax^), Akaike information criterion (AAIC), and 
Bayesian information criterion (ABIC) relative to ACDM 
value, and p-values estimated from the likelihood ratio test 
(LRT) statistic. 


Data 

ACDM 

Scaling-DE 

Spike-DE 

Lowlike 

2014.578 

2014.178 

2014.092 

Commander 

-7.304 

-7.471 

-8.096 

CamSpec 

7795.773 

7796.223 

7777.669 

Leasing 

9.892 

9.190 

9.881 

DR7 

0.858 

0.620 

0.439 

DR9 

0.431 

0.603 

0.812 

6dF 

0.019 

0.034 

0.036 

Wiggle Z 

0.047 

0.024 

0.021 

Total y^ 

9814.295 

9813.400 

9794.753 

Cl 

<1 


-0.895 

-19.542 

p-value (LRT) 


0.3441 

6.148 X lO”"' 

AAIC 


1.105 

-11.542 

ABIC 


6.982 

11.968 


over the other by the current observations based on the 
statistical criteria such as the likelihood ratio test and the 
Akaike and Bayesian information criteria that are widely 
used in the model selection. 

Table Hill lists the separate chi-square (y^) values for 
each likelihood data set used in the parameter estima¬ 
tion for the best-fit ACDM, scaling-DE, and spike-DE 
models (see Table U for the best-fit values). The negative 
chi-square values for CommEinder likelihood data appear 
due to the arbitrary normalization of log-likelihood in 
the CosmoMC software [1^. We note that the best-fit 
cosmological parameters in the presence of a dark en¬ 
ergy spike near the radiation-matter equality gives a far 
smaller chi-square value than those of the ACDM model 
by about five times the number of new free parameters of 
the spike-DE model (A, B, log^g ac, a) with a difference 
- Xacdm = -19.542, which is a significant 
improvement of data-fitting. 

The three dark energy models considered here are 
nested in the sense that the ACDM and scaling-DE mod¬ 
els are special cases of the spike-DE model. In this case, 
we can apply the likelihood ratio test (LRT) as a model 
selection method, where the null model is the ACDM 
model and the alternative model is the scaling-DE or 
spike-DE model UM- The test statistics is defined as 
the twice the natural logarithm of the ratio of likelihoods 
of the null and alternative hypotheses (models) and is 
equivalent to a difference of chi-square relative to the 
ACDM one, 


Q = 21n 


-^(.ffACDM|.P) 

CiH\D) 


Ax" 


(7) 


where C{H\D) indicates the maximum likelihood of the 
alternative model (H) given the data (D), and likewise 
for the null model (L^acdm)- The LRT statistic is a com¬ 
putationally cheap version of the Bayes factor which pro¬ 
vides a criterion for penalizing models with more param¬ 
eters based on the Bayesian theory 0- The test statistic 
Q can be approximated as the y^-distribution with de¬ 
grees of freedom {df) defined as the additional number 
of parameters of the nesting model {df = 4 for spike-DE 
model, and df = 1 for scaling-DE model). The p-value, 
the probability that the null hypothesis is supported by 
the observational data over the alternative one, is calcu¬ 
lated from the cumulative y^-distribution and presented 
in Table Ell We find that the p-value for the spike-DE 
model as alternative is quite small (p = 6.1 x lO”'^), sug¬ 
gesting a strong preference to the spike-DE model against 
the ACDM model. 

The Akaike information criterion (AIC) is defined as 


AIC = -2lnC + 2k, (8) 

where k is the number of free parameters of the model 
considered. If the alternative model gives a smaller AIC 
compared with the null (ACDM) model, it is ranked as 
a better model because the discrepancy with the true 
model is considered to be smaller. It is generally accepted 
that the AIC difference of 5 or more gives a strong evi¬ 
dence supporting the model with smaller AIC value (see 
Q for the reliability of the AIC method in cosmologi¬ 
cal model selection). The differences of AIC relative to 
the ACDM model (AAIC = Ay^ -I- 2df) are listed in 
Table ED The scaling-DE model has a positive value of 
AAIC = 1.1, which means that introducing the scaling 
dark energy without a spike into the ACDM model does 
not improve the fit much. On the other hand, the neg¬ 
ative value of AAIC = —11.5 for the spike-DE model 
suggests that the alternative model of dark energy with 
early scaling behavior and a dark energy spike (near the 
radiation-matter equality) is strongly favored by the cur¬ 
rent cosmological observations over the ACDM model. 

As an alternative to the AIC, the Bayesian information 
criterion (BIC) is often used for model selection, which 
assigns a conservative penalty for large sample size. The 
BIC is defined as sain 

BIC =-21n/:-kfclnA, (9) 

where N is the number of data points. We set N = 2637 
for Planck-l-BAO data sets (31 x 4 -|- 48 -f 2451 -I- 8 -I- 6 
for Lowlike [TT, TE, EE, BB], Commander, CamSpec, 
leasing, and BAO data, respectively) to calculate the 
difference of BIC relative to the ACDM value (ABIC = 
Ay^ + d/lnA), which are listed in Table Ell Contrary 
to the AIC result, the spike-DE model has a positive 
value of ABIC = 12.0. The strong evidence supporting 
the presence of dark energy spike by the AIC has disap¬ 
peared, since the BIC penalizes complex models by the 
large number of data points as in the CMB observation. 
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FIG. 3: Top: The CMB temperature angular power spec¬ 
tra of the best-fit ACDM (black), scaling-DE (grey), and 
spike-DE models (red curves). The angular power spec¬ 
trum is given as Di = liJL l)C'i/(27r) in unit. The 
Planck 2013 data points are shown with blue dots with er¬ 
ror bars. Small panels inside magnify regions of i ~ 800 
and 1800. Middle and Bottom: the difference between obser¬ 
vation and best-fit model prediction of band power divided 
by the measurement error (at) for each .^-bin (denoted as b), 
Zb — {Db^ohs—Db,inod)/o'b, and the sum of contribution due to 
Zb within the bin, Xb ~ Zb{A£)b. For comparison, the quan¬ 
tities Zb and Xb estimated from the Planck 2015 data and 
the best-fit ACDM model prediction (constrained with Plank 
2015 TT-|-LowP+Lesing data) are presented (green open cir¬ 
cles; 0]). 


In the context of BIG, the scaling-DE model with df = 1 
is preferred over the spike-DE model {df = 4). 

In summary, comparison of the maximum likelihoods 
of spike-DE and ACDM models according to LRT and 
AIC suggests that the spike-DE model is strongly pre¬ 
ferred over the ACDM one while the BIC still indicates 
the observational data supports the simple ACDM model 
over others. From the definition of AIC and BIC we see 
that the AIC is inclined to select the model that bet¬ 
ter fits to the data while the BIC selects a simpler model 
with less parameters. Apart from the model selection be¬ 
tween ACDM and spike-DE models, at least we conclude 
that the spike-DE model fits the observational data far 
better than ACDM model with the different cosmological 
parameter estimation. 

According to Table IIIIl the spike-DE model improves 
fitting to the CamSpec high i temperature likelihood data. 
Figure [3] verifies that the chi-square decrease is mainly 


due to the better fitting of Planck temperature power 
spectrum data around the third {£ « 800) and sixth 
{£ Ri 1800) acoustic peaks where strong residuals rela¬ 
tive to the best-fit ACDM model are seen. In the middle 
and bottom panels of Fig. [3l we plot the difference of 
the observed and the model-predicted band power spec¬ 
tra {Db) normalized with the measurement error (uf,) for 
each £-bin (here denoted as b), Zb = {Db,ohs— Db,mod)/o’b, 
and the sum of contribution due to Zf within the bin, 
= Z^{Al)b, which approximates the chi-square con¬ 
tribution for each bin. The model band power Db^mod is 
the average of the CMB angular power spectrum = 
£{£ -I- l)Q/(27r) predicted by a model within a specified 
bin. We use the same £-bins used in the Planck team’s 
analysis. It was originally reported that the strong resid¬ 
uals seen around the third {£ r; 800) and fifth {£ rs 1300- 
1500) acoustic peaks are real features of the primordial 
CMB sky [0. We note that the best-fit spike-DE model 
significantly alleviates the strong residuals around the 
third {£ r: 800) and sixth {£ Ri 1800) acoustic peaks ob¬ 
served in the best-fit ACDM model (red and black dots). 
However, the residual around the fifth peak {£ Ri 1300- 
1500) still remains in both models. 


V. DISCUSSION AND CONCLUSION 

In this paper, we investigate the observational effect of 
the early episodically dominating dark energy which ac¬ 
commodates a dark energy spike, a sudden transient vari¬ 
ation in dark energy density, together with early scaling 
behaviors and late-time acceleration. 

The dark energy model with a spike (spike-DE model) 
near the radiation-matter equality era improves the fit 
to the Planck CMB temperature power spectrum data 
around the third {i Ri 800) and sixth {£ r; 1800) acous¬ 
tic peaks. Comparing the likelihood distributions based 
on the maximum likelihood ratio test and the Akaike 
information criterion as the statistical model selection 
methods, we find that the spike-DE model is strongly 
favored by observations over the conventional ACDM 
model. Furthermore, the spike-DE model provides the 
different cosmological parameter estimation with tighter 
constraints on matter density and Hubble constant (see 
Fig. and Table HIl) . However, the strong evidence sup¬ 
porting the presence of the dark energy spike disappears 
based on the the Bayesian information criterion which 
assigns a conservative penalty to the model with a large 
number of parameters. 

We have checked that including high-£ CMB data ob¬ 
served by the South Pole Telescope and the Atacama 
Cosmology Telescope [H, [30 or excluding the tensor- 
type perturbation do not affect our main results. Besides, 
we infer that the foreground and instrumental calibration 
parameters do not play a major role in improving the fit 
to the data. If they do, a significant reduction of chi- 
square should be seen in the case of ACDM model, too. 

Very recently, the Planck 2015 data have been pub- 



























licly available. The main scientific conclusions of Planck 
2015 data analysis are consistent with the previous re¬ 
sults, with cosmological parameters deviating less than 
0.7(7 Q. As shown in Fig. [31 the strong residuals around 
£ ~ 800 and ^ ~ 1800 in the Planck 2013 temperature 
power spectrum data are not observed in the Planck 2015 
data; the deviation from the best-fit ACDM model pre¬ 
diction becomes quite smaller (see green open circles in 
the middle and bottom panels of Fig. O. Therefore, the 
success of the spike-DE model in improving the fit to 
the Planck 2013 temperature power spectrum at those 
regions is not expected any more in the recent Planck 
2015 data. However, the presence of strong residuals at 
£ = 400-500 and £ ~ 1200 in the Planck 2015 data still 
leaves open the possibility that the new data are fitted by 
another candidate model of dark energy far better than 
ACDM model. 


Through an example of the dark energy spike model, 
we emphasize that the alternative cosmological parame¬ 
ter estimation is allowed in the Einstein’s gravity, with 
even better fitting of the same observational data than 
the conventional ACDM model. 
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